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Abstract 



Thepp — * nK + Y> + reaction is a very good isospin 3/2 filter for studying A ++ * decaying 
to K + T, + . With an effective Lagrangian approach, contributions from a previous ignored 
sub--fT + £ + -threshold resonance A ++ *(1620)l/2~ are fully included in addition to those 
already considered in previous calculations. It is found that the A ++ *(1620)l/2~ reso- 
nance gives an overwhelmingly dominant contribution for energies very close to threshold, 
with a very important contribution from the t-channel p exchange. This solves the prob- 
lem that all previous calculations seriously underestimate the near-threshold cross section 
^ ' by order(s) of magnitude. Many important implications of the results are discussed. 



PACS: 13.75. Cs.; 14.20.Gk.; 13.30.Eg. 



1 Introduction 

The spectrum of isospin 3/2 A ++ * resonances is of special interest since it is the most exper- 
imentally accessible system composed of 3 identical valence quarks. However, our knowledge 
on these resonances mainly comes from old ttN experiments and is still very poor [1, 2, 3]. 
Many model-predicted A* resonances have not been established or ever observed [3, 4]. A 
possible reason is that these "missing" baryon resonances have too weak coupling to ttN to be 
observed in the nN experiments. Searching for these "missing" A* resonances from other pro- 
duction processes is necessary. A possible new excellent source for studying A ++ * resonances 
is pp — > nK + Y, + reaction, which has a special advantage for absence of complication caused by 
iV* contribution because of the isospin and charge conversation. 

At present, little is known about the pp — > nK + H + reaction. Experimentally there 
are only a few data points about its total cross section versus energy [5, 6]. Theoretically a 
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resonance model with an effective intermediate A ++ *(1920) resonance [7, 8, 9] and the Jiilich 
meson exchange model [10] reproduce the old data at higher beam energy [5] quite well, but 
their predictions for the cross sections close to threshold fail by order of magnitude compared 
with very recent COSY-11 measurement [6]. Here we restudy this reaction by examining various 
possible sources for the very strong near-threshold enhancement. 

In the resonance model with an effective intermediate A ++ *(1920) resonance [7, 8, 9], 
the 3/2 + A ++ *(1920) resonance decays to in relative P-wave and is suppressed at lower 

energies. To reproduce the near-threshold enhancement for the pp — > nK + T l + reaction, a 
natural source could be some l/2~ A ++ * resonance(s) at lower energy decaying to K + H + in 
relative S-wave. In the Jiilich meson exchange model [10] , only 7r and K exchanges are included 
with ir + p — ► K + H + and KN — > KN amplitudes taken from the relevant existing data. Its 
failure at lower energy indicates that other meson exchange might be important. Following 
the logic, we find a natural source for the near-threshold enhancement of the pp — > nK + H + 
reaction coming from p + p -> A++*(1620)(l/2-) -> The well established A++*(1620) 

has an unusual large decay branching ratio to the pp + [1, 11] and its coupling to K + T, + can be 
obtained from its coupling to 7r + p with SU(3) symmetry. Possible n-S + final state interaction 
are also studied. 

In next section, we will give the formalism and ingredients for our calculation. Then 
numerical results and discussion are given in Sect. 3. 



2 Formalism and ingredients 



The basic Feynman diagrams for the pp — > nK + Y, + reaction are depicted in Fig. 1. Besides the 
A ++ *(1920)3/2 + considered in the previous resonance model [7, 8, 9], A ++ *(1620)l/2~ reso- 
nance is added in our calculation. Both tt + and p + exchanges are considered for the production 
of A ++ * resonances. A Lorentz covariant orbital-spin(L-S) scheme [12] are used for the effective 
A*ttN, A*pN, and A*KE vertices. 
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(a) (b) 
Figure 1: Feynman diagrams for pp — > nK + Y> + reaction. 
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2.1 Meson-Baryon-Baryon (Resonances) vertices 

The relevant effective Lagrangian densities for the meson-NN vertices are the standard ones [7] : 

ArVTV = -WitNnN^T ■ 7rN, (1) 

C p nn = -9 p nnN(^ + -H-a^d")? ■ ?N. (2) 

Zttin 

And the relevant off-shell form factors for the 7r + -NN and p + -NN vertices are taken as the same 
as the well known Bonn potential model [13]: 



A 2 2 

pNN( n 2 \ _ ( A M ~ m M \n (o\ 
1V M <iM 



with n=l for 7T + -meson and n=2 for p + -meson. qM, tum and Km are the 4-momentum, mass 
and cut-off parameter of the exchanged- meson (M), respectively. Parameters are taken as 
commonly used ones [13, 7, 8]: gl NN /A-K = 14.4, = 1.3 GeV, g 2 pNN /A-K = 0.9, A p = 1.85 
GeV, and k — 6.1. 

To calculate the amplitudes for diagrams in Fig. 1 with resonance model, we also need to 
know interaction vertices involving A* resonances. In Ref. [12], a Lorentz covariant orbital-spin 
scheme for N*NM couplings has been described in detail. The scheme can be easily extended 
to describe the A*ttN, A*pN and A*i^S couplings that appears in the Feynman diagrams 
showing in Fig. 1. 

Firstly, for A ++ * — > 7r + p, it is well known that there is only one possible L-S coupling 
for the 7r + p final state of each A ++ * decay. Since the nucleon has spin-parity | + and pion 

has spin-parity 0~, with the parity and angular momentum conversation, A ++ *(1620)(| ) and 
A ++ *(1920)(| + ) can only decay to 7r + p in S-wave and in P-wave, respectively. The effective 
A ++ *7r + p couplings in the covariant orbital-spin scheme can be easily obtained as : 

A ++ *(1620)(i )^7T + p : 9 A* (1620) NnUpU*, (4) 

3 + 

A++*(1920)(- ) ^n + p : g a* (1920) NnU p u^ : (5) 

with Up, and as the Rarita-Schwinger spin wave functions for the proton, A ++ *(1620) 
and A ++ *(1920), respectively. 

Secondly, for A ++ * — > K + H + , the effective A ++ *K + H + couplings can be easily obtained 
with SU(3) symmetry by simply replacing n + and p with K + and E + , respectively. 

Thirdly, for A ++ * — > p + p, unlike pion with spin 0, here p has spin 1. For a A* with spin 
I there are two independent L-S couplings conserving parity and total angular momentum; for 
a A* with spin larger than |, there are three independent L-S couplings [12]. Here, in our 
calculation, we only consider the A ++ *(1620)(| ) resonance which was found to have large 
decay branch ratio to pN in relative S-wave [1, 11]. Following Ref. [12], the effective coupling 
is obtained as 



A++*(1620)(- ) - p + p : ^1900)^75(7^ - ^f^)u^, (6) 
with p* and the four momentum of the A* and the polarization vector of p meson. 
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For the relevant vertices involving A*, since only S- and P- waves are involved, we use the 
monopole form for the off-shell form factors : 



F M ( q ) = % m K (7) 

with A; = A; = 1.3 GeV. 

2.2 Coupling constants for A ++ * resonances 

The relevant A*-baryon-meson coupling constants are determined either from experimentally 
observed partial decay widths or SU(3) symmetry. 

The general formula for the partial decay width of a A* resonance decaying to a nucleon 
and a pion is as the following 

d 3 pN tti n d 3 p n 1 



where M.A*^m represents the total amplitude of the A* decay to a nucleon and a pion, the 
Pa*, Pn and p n are the four momentum of the three particles, En and E n are the corresponding 
energies. 

For A*(1620)(i") -> Ntt, the partial decay width can be calculated with the amplitude 
given by Eq. (4) to be 

„ _ 9l* (im N,( m N + E N )pT 

1 A*(1620)->iV7r - T—rZ (9 J 

47TMa*(1620) 

with 



P C N = 



( M A*(1620) ~ ( m N + "V) 2 )(^A*(1620) ~ ^ ~ m ^ nm 

\ AM 2 ' ^ ' 

\ ^ Ki A*(1620) 

En = s]{pT) 2 + ™n- (11) 

For A* (1920) (| ) -> Ntt, the partial decay width can be calculated with the amplitude 
given by Eq.(5) to be 

<?A*(1920)^Kv + EN^pTf , 10 . 

1 A*(1920)A% - — M \M) 

1^7TMa*(1920) 

with p^ here the momentum of the nucleon in the rest frame of A* (1920). 

For A* — > KE, the formulae are basically identical as for A* — > nN with the replacement 
of it to kaon and proton to E. 

For A*(1620)(|") -> Np -> Nttit , the partial decay width of A*(1620) -> Np in relative 
S-wave via p decay into two pions is given by 

r T~n -d 3 p 1 m l d 3 p 2 1 d 3 p 3 1 ^4/,^ \ 

r A ni620)^ = \M A ^ Np ^\ J^y^J^ W 2 j2^T 2Es { } (Ma * (1620) ~ Pl -P*-P 3 > 

(13) 
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Table 1: Relevant A* paramters. 



Resonances 


Width 


Decay 


Branching 


0747T 




(MeV) 


Channel 


ratios (%) 




A*(1620) 


150 


TiN 


25 


0.06 






pN 


14 


0.37 






KH 




0.06 


A*(1920) 


200 


ttN 
KY, 


12.5 
2.1 


0.18(GeV- 2 ) 
0.12(GeV- 2 ) 



where M.£s^,n p -+Ni™ represents the total amplitude of A*(1620)(| ) — > Np — > Ntttt. p ± , m 1; 
and Ei stand for the four momentum, mass, and energy of the nucleon; p2, P3, and E2, E% 
stand for the four momentum and energy of the final two pions, respectively. In the amplitude 
calculation, the amplitude for p — > tytt is taken as M. fm -w = g P n^p,(p, m p){P2 ~ P3) with e p the 
polarization vector of p meson and g 2 /Air = 2.91 [14]. 

With the experimental branching ratios [1] and above formulae, we can obtain all relevant 
coupling constants as summarized in Table 1. Because A ++ *(1620)(| ) has its mass below 
K + T, + threshold, we take 5 , A*(i62o)ftrs/^ 7r = 5 , A*(i62o)A r 7r/^ 7r = from SU(3) relation. 

2.3 Propagators 

In our amplitude calculation of Feynman diagrams in Fig. 1, we also need propagators for the 
pion, p + meson, and the intermediate A ++ * resonances with half-integer spin. For the pion and 
the p + meson, the propagators are : 

GM = -r^-ii ( 14 ) 

J2 



G7(? P ) = , 3 ) ( 15 ) 



q 2 p -m 2 p 

with q n and q p the four momenta of pion and p + meson, respectively. m n and m p are the 
corresponding masses. 

The propagators for the A ++ * resonances with half-integer spin can be constructed with 
their projection operators and the corresponding Breit-Wigner factor [15]. For the spin-(n+|), 
the propagator can be written as 

where 1/(-Pr — M R + iMrYr) is the standard Breit-Wigner factor; Mr, Pr and Tr are the 
mass, four momentum and full width of the resonance, respectively. p( n+ 2) is the projection 
operator. Explicitly, 

1 -f-p + M R 
P2{P)= 2M R ' (17) 



2M, 



R 



R ft 



(18) 
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for A*(1620)(i ) and A*(1920)(§ + ), respectively. 



2.4 Amplitude and total cross section for pp — > nK + T l + 

The full amplitude in our calculation for the pp — > nK + H + reaction is composed of three parts 
corresponding to A ++ *(1620) production from 7r + exchange, A ++ *(1620) production from p + 
exchange and A ++ *(1920) production from 7r + exchange, respectively. Explicitly, 

M = 7W(A++*(1620), tt + ) + 7W(A ++ *(1620), p + ) + 7W(A ++ *(1920), tt + ). (19) 

Each amplitude can be obtained straightforwardly with effective couplings and propagators 
given in previous sections by following the Feynman rules. Here we give explicitly the amplitude 
_M(A++*(1620),7r + ), as an exam pi e) 

.M(A ++ *(1620),7T + ) = Ux(px,Sx)g A * Y ;KG*(p*)gA*N7rU N (p 1 ,Si)G 7r (q 7r ) x 

u N (p n , s n )V2g nNN -f 5 u N (p 2 , s 2 ) 

+ (exchange term with p 1 <-> p 2 ) (20) 

with se, s n , si, s 2 the spin projection of S + , neutron in the final state and two initial protons, 
respectively. G*(p*) and G w are the propagators for A*(1620) resonance and the exchanges 
pion. un and uy, are the Dirac wave functions of the nucleon and the S. p\ and p 2 represent 
the four momenta of the initial protons. 

Then the calculation of the invariant amplitude square \M. | 2 and the cross section a{pp — > 
nK + Y, + ) is straightforward. 



a ( lm P ST^ \ KA \2 m nd i p n d i p K mYd 6 pY 4 

da(pp ^ nK + X + ) = --jf ^ 1^1 — g ^ ^ 6 (Pi + ^ - - - Ps) (21) 

si,S2,s n ,s B n A" E 

with the flux factor 



F=(27r) 5 v /(p 1 -p 2 ) 2 -m4. (22) 

Since the relative phases between three parts in the amplitude of Eq.(19) are not known, the 
interference terms between different parts are ignored in our concrete calculations. 



2.5 Final State Interaction 

For three particles in the final state, the interaction between K + and S + is dominated by 
the s-channel A ++ * resonances and is included in our resonance isobar model calculation; the 
S-wave interaction between K + and nucleon is known to be weak and repulsive [16] and does 
not play significant role [17, 18]. So we only need to consider n-£ + final state interaction. 

In Ref. [6], the total cross-section was measured at the COSY-11 spectrometer at excess 
energies Q=13 MeV and Q=60 MeV. At such near-threshold energies, the n-S + final state 
interaction could play very significant role according to the experience in studying a similar 
process pp — > pK + A both experimentally [17] and theoretically [18]. To study possible influence 
from the n-S + final state interaction, we include it in our calculation with the Watson-Migdal 
approach [19] by factorizing the reaction amplitude as: 

A = M( PP -> nX + E+)T nS (23) 
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where M.(pp — > nK + H + ) is the primary production amplitude as discussed above, T n s describes 
the n-S + final state interaction, which goes to unity in the limit of no FSI. For the near-threshold 
FSI, the enhancement factor T„ s can be expressed similarly as for p-A FSI in the study 
of pp — > pK + A process [18]: 

T nS = »±£ (24) 
q — la 

where q is the internal momentum of subsystem, and the a and (3 are related to the 

effective-range parameters via 

a/3 a + 

Lacking the knowledge on the n-S + interaction, in the present work we assume similar a 
and P values for as for the p-A [18], i.e., 

a = -75.0 MeV, (3 = 200.0 MeV. (26) 



3 Numerical results and discussion 

With the formalism and ingredients given above, the total cross section versus the kinetic energy 
of the proton beam (Tp) for the pp — > nK + T, + reaction is calculated by using a Monte Carlo 
multi-particle phase space integration program. The results for Tp between 1.5 and 6.0 GeV 
are shown in Fig. 2 together with experimental data [5, 6] for comparison. 




Figure 2: Total cross section vs Tp for the pp — > nK + H + reaction from present calculation (solid 
curves) compared with data [5, 6]. Left: dot-dashed, dashed and dotted curves for contributions 
from A*(1620)(7r + exchange), A*(1620)(p + exchange) and A*(1920)(7r + exchange), respectively. 
Right: contribution of effective A* (1920) scaled by a factor 5 with FSI (dotted curve) and 
without FSI (dot-dashed); Full calculation with FSI (solid) and without FSI (dashed). 
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In the left figure of Fig. 2, contributions from A*(1620)(7r + exchange), A*(1620)(p + ex- 
change) and A*(1920)(7r + exchange) are shown separately by dot-dashed, dashed and dotted 
curves, respectively. The contribution of the A* (1620) production by the p + exchange is found 
to be very important for the whole energy range. The contribution of the A* (1620) resonance 
is found to be overwhelmingly dominant over A* (1920) contribution for the two lowest data 
points close to the threshold. This gives a natural source for the serious underestimation of the 
near-threshold cross sections by previous calculations [7, 8, 9, 10], which have neglected either 
A* (1620) resonance contribution [7, 8, 9] or p + exchange contribution [10]. The solid curve in 
the figure is the incoherent sum of the three contributions and reproduces the data quite well 
for the whole energy range. The remained slightly underestimation for the two lowest data 
points may be improved by the interference term between p + exchange and ir + exchange for 
the A*(1620) production. 

To show the effect from the n-S + FSI, we give the result without including the FSI factor 
by the dashed curve in the right figure of Fig. 2. Comparing dashed curve with the solid curve 
which includes the FSI factor, we find that the FSI enhances the total cross section by a factor 
of about 3 for the two lowest data points. So the FSI is indeed making a significant effect at 
energies close to threshold. But it does not change the basic shape of the curve very much. 
In previous calculations [7, 8, 9], only A* (1920) contribution are considered with a free scaling 
parameter to fit the data. In Fig. 2 (right), we also show the results from only A*(1920)(7r + 
exchange) scaled by a factor 5 for comparison. It reproduces the data for Tp above 2.8 GeV 
quite well, but underestimates the two lowest data points by orders of magnitude no matter 
whether including the FSI (dotted curve) or not (dot-dashed curve). 

In Fig. 3, we give our model prediction of the Dalitz Plot and all relevant mass spectra 
for the reaction at T P = 2.8 GeV which can be achieved by the proton beams at Lanzhou 
Cooler Storage Ring (CSR) and COSY at Julich. A strong K + H + near-threshold enhancement 
is predicted together with a relative weak broad peak around 1920 MeV in the K + H + invariant 
mass spectrum and some nS + near-threshold enhancement. The prediction may be checked by 
future experiments at Lanzhou CSR with the scheduled 4n hadron detector [20] or at COSY 
with the newly installed WASA-at-COSY detector [21]. 

The result on the dominant role of the A ++ *(1620)l/2~ resonance for the near-threshold 
cross section of pp — > nK + Y> + reaction has many important implications. 

(1) It gives another example that sub-threshold resonances can make extremely important 
contributions and should not be simply ignored. Many calculations were used to consider 
only the resonances above threshold, such as previous calculations [7, 8, 9] for pp — > pK + A 
and pp — > nK + H + . In the case of pp — > pK + A, the sub--fT + A-threshold resonance A r *(1535) 
is proposed to play a significant role for the near-threshold energies [22]. Comparing with 
recent experimental measured Dalitz plot by COSY-TOF Collaboration, a model fitting without 
including the contribution of the iV*(1535) visibly underestimates the part near KA threshold 
[23]. The important role of the iV*(1535) for the KA production is most clearly demonstrated 
by the BES data on J/ip — > pK + A [22, 24]. There are several more examples from J/i[) 
decays showing the importance of contribution from sub-threshold particles, such sub-7riV- 
threshold nucleon pole contribution in J/ip —> pmr + [25, 26], sub-i^^-threshold contribution 
in J/i/j — > KKn and sub-cj7r-threshold contribution in J/ip — » unn [27]. 

(2) Since the pp — > nK + T l + and pp — » pK + A reactions are the basic inputs for the 
strangeness production from heavy ion collisions [28, 29], the inclusion of the sub-threshold 
A* ++ (1620) and iV*(1535) contributions may be essential for such studies. 

(3) The t-channel p exchange may play important role for many meson production pro- 
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cesses in proton-proton collisions and should not be ignored. 

(4) The extra-ordinary large coupling of the A* (1620) to pN obtained from the n + p — > 
Nnn [1, 11] seems confirmed by our present study of the strong near-threshold enhancement 
of pp — > nK + T l + cross section. Does the A* (1620) contain a large pN molecular component 
or relate to some pN dynamical generated state ? If so, where to search for its SU(3) part- 
ners ? Sarkar et al. [30] have studied baryonic resonances from baryon decuplet and psudoscalar 
meson octet interaction. It would be of interests to study baryonic resonances from baryon 
octet and vector meson octet interaction. In fact, from PDG compilation [1] of baryon reso- 
nances, there are already some indications for a vector-meson-baryon SU(3) decuplet. While 
the A* (1620) l/2~ is about 85 MeV below the Np threshold, there is a £*( 1750) 1/2" about 
70 MeV below the NK* threshold and there is a S*(1950)? ? about 60 MeV below the AK* 
threshold. If these resonances are indeed the members of the l/2~ SU(3) decuplet vector- 
meson-baryon S-wave states, we would expect also a f2*l/2~ resonance around 2160 MeV. All 
these baryon resonances can be searched for in high statistic data on relevant channels from 
vector charmonium decays by upcoming BES3 experiments in near future. 
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